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G> ■ Abstract 

We perform a global analysis of the lepton-quark neutral current data on the low scale gravity model, 
which arises from the extra dimensions. The global data include HERA neutral current deep-inelastic 
scattering, Drell-yan production at the Tevatron, and fermion pair production at LEPII. The Drell-yan 
production, due to the large invariant mass data, provides the strongest constraint. Combining all data, 
the effective Planck scale must be larger than about 1.12 TeV for n = 3 and 0.94 TeV for n = 4 at 
95%CL. 
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' 1. Introduction 

o ■ 

Q>^ i A solution to the gauge hierarchy, based on extra spacial dimensions, was recently proposed by Arkani- 

f~i ' Hamed, Dimopoulos and Dvali W. They assumed the space-time is 4 + n dimensional, with the standard 

I ' model (SM) particles living on a brane. While the electromagnetic, strong, and weak forces are confined 

O-t. . . 

, to this brane, gravity can propagate in the extra dimensions. To solve the gauge hierarchy problem they 

proposed the "new" Planck scale Ms is of the order of TcV with very large extra dimensions. The size R 
■ of these extra dimensions can be as large as 1 mm, which corresponds to a compactification scale R^^ as 

5— ( ' low as 10"'* eV. The usual Planck scale Mq — ^/VGn 1-22 x 10*^ GeV is related to this effective Planck 

scale AIs using the Gauss's law: 

R"M^+^r^M^, (1) 

where n is the number of extra dimensions. For n = 1 it gives a large value for R, which is already ruled out 
by gravitational experiments. On the other hand, n = 2 gives R<1 mm, which is in the margin beyond the 
reach of present gravitational experiments. 

The graviton including its excitations in the extra dimensions can couple to the SM particles on the brane 
with an effective strength of 1/Ms (instead of I/Mq) after summing the effect of all excitations collectively, 
and thus the gravitation interaction becomes comparable in strength to weak interaction at TeV scale. Hence, 
it can give rise to a number of phenomenological activities testable at existing and future colliders So 
far, studies show that there are two categories of signals: direct and indirect. The indirect signal refers to 
exchanges of gravitons in the intermediate states, while direct refers to production or associated production 
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of gravitons in the final state. Indirect signals include fermion pair, gauge boson pair production, correction 
to precision variables, etc. There are also other astrophysical and cosmological signatures and constraints 
1^ . Among the constraints the cooling of supernovae by radiating gravitons places the strongest limit on the 
effective Planck scale Ms of order 50 TeV for n = 2, which renders collider signatures for n = 2 uninteresting. 
Thus, we concentrate on n > 3. 

In this work, we perform a global analysis of the lepton-quark neutral current data on the low scale gravity 
model. The global data include HERA neutral current deep-inelastic scattering, Drell-yan production at the 
Tevatron, and total hadronic, bb and cc pair cross sections at LEPII. In addition, we also include the leptonic 
cross sections e+e^ /i^ fi^ ,t^t^ at LEPII in our analysis. The t^-N scattering data from CCFR and 
NuTeV have been shown by Rizzo |^ to be very insignificant in constraining the low scale gravity, and so we 
shall not include this data set in our analysis. We shall see that the Drell-yan production, due to the large 
invariant mass data, provides the strongest constraint among the global data. By combining all data, the 
effective Planck scale Mg must be larger than about 1.12 TeV for n = 3 and 0.94 TeV for n = 4 at 95%CL. 
This is our main result. 

The organization of the paper is as follows. In the next section, we shall describe each set of data used 
in our global analysis and derive the effect of the low scale gravity. In Sec. Ill, we give the numerical results 
and interpretations for our fits, from which we can draw the limits on the effective Planck scale. 



2. Global Data 



Before we come to each data set, let us first give a general expression for the square of amplitude for 

e^e+ qq: 

2 



= (IMllP + iM^flp) + 4<2 (iMiflp + \Mrl\ 



2n' 



S / 



Sne^QeQ: 



T (u~tf 



« Ml 



87re^ 



(2) 



sin^ 6'w cos2 6'w s - Ml 
where s,t,u are the usual Mandelstam variables, e = \/47ra is the electromagnetic coupling constant, Qj 
is the electric charge of the fermion /, 0^ is the Weinberg mixing angle, g[ and gl are, respectively, the 
axial-vector and the vector Z couplings to the fermion /. The reduced amplitudes Maf3 {a,f3 = L,R) are 



Ma/3 = 



6 QeQq 



sin^ 6'w cos2 6'w s-Mj. 



(3) 



where 



gi = Tsf - Qf sin^ 6'w , 9r^ ^Qf sin^ ^'w 
9l^{9l+9i)/2, gl^{gl-g{,}/2 
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In the derivation, we have foUowed the notation in Han et al. ^ and we have taken Mg ^ s, |u| and in 
this case the propagator factor D{s) = D{\t\) = D{\u\), which is given by 



= ^ x^, (4) 



where the factor J- is given by 



{log ( for n — 2 , 
^ V ^ ^ (5) 
^ for n > 2 . 

For n > 3 ^ is always positive. The amphtude square for the crossed channels 9 — > Q can be 
obtained using the crossing symmetry. 

2. a HERA neutral-current data 

Both HI and ZEUS have measured the neutral-current deep-inelastic scattering cross sections at high-Q^ 
region. Despite the excess in cross section reported in early 1997, the 1997 data alone agreed well with the 
SM. The double differential cross section for e+p e^X is given by 



dydQ 




sm^ t/w cos 



^2 ^ fQ^ 



+ 2^^(^)i7 [mI) [f) + (6) 

where = sxy is the square of the momentum-transfer, ^/s — 300 GeV, and fq/g{x) are the parton 
distribution functions. The reduced amplitudes Maf3 are given in Eq. (|) with s replaced by -Q^. The 
last term in the above equation comes from the additional channel e^g — > e^g allowed by the low scale 
gravity interactions. The channels X^n^ 9 ^ ^ 9 

also included. The data of ZEUS [|, including 
systematics, are recently published while the HI data used in our analysis are shown in Table |l|. 

2.b Drell-yan at the Tevatron 

CDF measured the double differential cross section (Pa/dMudy (£ = e,/i) for the Drell-yan production, 
where AIu and y are, respectively, the invariant mass and the rapidity of the lepton pair. Essentially, CDF 
measured the cross section in the range — 1 < y < 1 and then average over y to obtain the double differential 
cross section. The double differential cross section for pp t^i" is given by 

Ml 



- Ay. 

dMudy 1 ^ 



727rs 



/g(a;i)/g(x2) (iMLip + lAfiflp + iMj^ip + \Mrr 



3 



(7) 



where xi,2 — Mue^^ / \fs, \fs = 1800 GcV, fq_g{xi) are the parton distribution functions, and the sum over 
ah possible qq pairs is understood. Note that after integrating the scattering angle cos 6** over the whole 
range in the center-of-mass frame, the interference terms proportional to [J- /Mg) vanish. The reduced 
amplitudes Map (a, P — L,R) are given in Eq. (^) with s replaced hy s — The QCD K factor is given 
hy K = I + ^^f^|(l + ^f-)- The second term of the above equation comes from qq ^ G ^ I'^i^ and the 
third one comes from gg G ^ i'^l". This K factor is, in principle, not valid for the gg channel, but it 
will not affect our calculation because we normalize our SM calculation to the expected number of events in 
each bin given by CDF when we deal with the CDF data (similarly, we normalize to the first bin of the DO 
data when we deal with DO data.) On the other hand, DO measured the differential cross section da/dMee 
with the rapidity range integrated. The Drell-yan data |^ we used in our analysis are shown in Table ^j. 

2.C LEPII fermion pair Cross sections 

LEPII has measured the leptonic cross sections, hadronic cross sections, and heavy flavor (6 and c) production. 
Since the low scale gravity interactions are naturally flavor-blind, we shall include the effects on all flavors 
here. For a massless q the expression for a{e^e~ — > qq) is given by 

a{e+e- ^ qq) - i^jj^ (|MllP + \Mlh\' + {MulI' + IMurI') + ^ (^j^)'} ' («) 

where the reduced amplitudes Map {a, (3 = L,R) are given in Eq. (||), and the QCD K factor is if = 
1 + ag/n + 1.409(as/7r)^ — 12.77 {as / n)'^ ■ The leptonic cross sections CT(e+e^ fj.'^ pi^ ,t^t^) do not have 
this K factor. The LEPII data |^ used in our analysis are given in Tables I and |. 

2.d Z/-N Scattering 

The recent measurements by CCFR and NuTeV Q are the most precise ones on the neutrino-quark couplings. 
Since the gravitons also couple to neutrinos, the measurements should, in principle, place a constraint on 
the scale Ms- However, the analysis by Rizzo in Ref. showed that the constraint coming from these low 
energy z/-N scattering is very weak. This is because the center-of-mass energies of these collisions are only of 
order of tens of GeV's, even though the neutrino beam energy is as high as a few hundred GeV. In addition, 
the effective operators induced by the low scale gravity interactions are at least of dimension eight. This is 
in contrast to the traditional four-fermion contact interactions, which are only of dimension six. Therefore, 
at such low center-of-mass energies in these collisions, the effect of low scale gravity is extremely minimal. 
We are not going to include these data in our global fit. 
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3. Fits 



For the fitting we follow the procedures in Refs. where the four-fermion contact interactions are analyzed 
with a similar but larger global data set. The difference is that Refs. Q include also the data sets from 
the low energy e-N scattering and the atomic parity violation, which constrain parity violating interactions. 
But the gravity is certainly parity conserving and thus it receives no restriction from these parity-violating 
experiments. 

The interference effect between the SM amplitude and the low scale gravity scales as T /Mg, while the 
pure low scale gravity scales as (JF/Af|)^. To linearize the fitting we use the parameter rj = T jMg. We 
use MINUIT to minimize the x^- The best estimates of r\ for each individual data set and the combined set 
are shown in Table |[ The SM fit gives a X — 94.10 for 118 d.o.f., while the fit with non-zero r\ gives a 

— 94.05 for 117 d.o.f. Since the fits do not show any evidence for new physics, we can then place limits 
on the r\ = T /Mg. The physical region of ?/ is 77 > 0, and so we define the 95%CL upper limit for 77 as 

95 = do Ul L (q] 

Jo T^iv) dr] 

where 



Viv) = exp ( - ) . (10) 



From ry+ we find the limits as Ms/T^^'^ = ?7+^^''. The limits on MsjT^^^ for each data set and the combined 
are also shown in Table ^. The combined limit is MsjT^I'^ > 939 GeV, or 

f 1120 GeV for n = 3 
Ms>{ (11) 
t 939 GeV for n = 4 

at 95%CL. 

It is the Drell-yan production at the Tevatron that gives the strongest constraint among all data. This 
is easy to understand because the s of the Drell-yan process is the largest of all. The new HERA data are 
now consistent with the SM and only cause a very slight pull of the fit, about one third of a from zero: 
see Table |. Both the "ALL w/o HERA" and the "ALL" fits agree well with the SM. We show in Fig. |l| 
the Drell-yan production for the SM and for the low scale gravity model with il/5 and n given in Eq. (p^), 
together with the DO data. 

To conclude, we have used the global lepton-quark neutral current data plus the leptonic cross sections 
at LEPH to constrain the low scale gravity interactions arising from the extra dimensions. The limit that 
we place on the effective Planck scale Ms is about 1.12 TeV for n = 3 and 0.94 TeV for n = 4 at 95%CL. 
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Figure 1: The differential distribution da/dM^e for Drell-yan production at the Tevatron for the low scale 
gravity model and for the SM. The data shown arc from DO. 
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Table 1: The measured iVobs and expected Ne^p number of events as a function of Q^i^ for HI at HERA. 



HI (£ = 37.04 pb"^) 



QLn(GeV^) 




N 


2500 


1297 


1276±98 


5000 


322 


336±29.6 


10000 


51 


55.0±6.42 


15000 


22 


14.8±2.13 


20000 


10 


4.39±0.73 


25000 


6 


1.58±0.29 
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Table 2: The Drell-Yan data by CDF and DO A'obs is the observed and TVoxp is the expected number of events. 



CDF 






e 






Mu (GeV) 


Nohs 




Nohs 


^ ''cxp 


50-150 


2581 


2581 


2533 


2533 


150-200 


8 


10.8 


9 


9.7 


200-250 


5 


3.5 


4 


3.2 


250-300 


2 


1.4 


2 


1.3 


300-400 


1 


0.97 


1 


0.94 


400-500 


1 


0.25 





0.27 


500-600 





0.069 





0.087 



DO 






Mu (GeV) 


a (pb) 


120-160 


1.93 


160-200 


0.49 1^:1^ 


200-240 


0.28 


240-290 


n l^ar +° 052 
0.066 -0.058 


290-340 


0.033 -0.030 


340-400 


+0.042 

0.057 -0.047 


400-500 


< 0.063 (0.039) 


500-600 


< 0.060 (0.037) 


600-1000 


< 0.058 (0.035) 



9 



Table 3: The LEPII hadronic cross sections and the ratio Ri, and Re- Note that the cross sections given by ALEPH 
have an angular cut of | cosd\ < 0.95. The two sets of DELPHI data for Rb,Rc are given at ^/s = 130 — 136 and 
161 - 172 GeV, respectively, crf^ and Rf^ are the SM expectations. 



(GeV) 


O-had (Pb) 


^h^d (pb) 


Rb 




Rc 




ALEPH 


130 


71.6 ± 3.96 


70.7 


0.176 ± 0.044 


0.190 






136 


58.8 ± 3.61 


57.3 


0.214 ± 0.050 


0.186 






161 


29.94± 1.84 


30.7 


0.159 ± 0.035 


0.175 






172 


26.4 ± 1.75 


25.1 


0.134 ± 0.037 


0.173 






183 


21.71± 0.74 


21.1 


0.176 ± 0.020 


0.171 






189 


19.59 ± 0.46 


19.3 










DELPHI 


130.2 


82.1 ± 5.2 


83.1 


0.174 ± 0.028 


0.182 


0.199±0.073 


0.225 


136.2 


65.1 ± 4.7 


67.0 










161.3 


40.9 ± 2.1 


34.8 


0.142 ± 0.024 


0.165 


0.314±0.055 


0.250 


172.1 


30.3 ± 1.9 


28.9 










L3 


130.3 


81.8 ± 6.4 


78 










136.3 


70.5 ± 6.2 


63 










140.2 


67 ± 47 


56 










161.3 


37.3 ± 2.2 


34.9 










170.3 


39.5 ± 7.5 


29.8 










172.3 


28.2 ± 2.2 


28.9 










OPAL 


130.25 


64.3 ± 5.1 


77.6 










133.17 






0.195 ± 0.041 


0.182 






136.22 


63.8 ± 5.2 


62.9 










161.34 


35.5 ± 2.2 


33.7 


0.162 ± 0.041 


0.169 






172.12 


27.0 ± 1.9 


27.6 


0.131 ± 0.047 


0.165 






183 


23.7 ± 0.81 


24.3 


0.190 ± 0.022 


0.163 






189 


22.1 ± 0.64 


22.3 


0.167 ± 0.014 


0.162 
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Table 4: The LEPII Icptonic cross sections for e+e M^^f , t^t ■ Note that the cross sections given by ALEPH 
have an angular cut of | cos^| < 0.95. a^^'^ are the SM expectations. 



(GeV) 




4f (pb) 


O-rr (pb) 


(pb) 


ALEPH 


130 


7.9 ± 1.2 


7.0 


10.9 ± 1.8 


7.3 


136 


6.9 ± 1.1 


6.1 


5.6 ± 1.3 


6.3 


161 


4.49± 0.70 


3.88 


5.75 ± 0.99 


4.01 


172 


2.64± 0.53 


3.32 


3.26 ± 0.75 


3.43 


183 


2.98± 0.25 


2.89 


2.90 ± 0.30 


2.98 


189 


2.66± 0.14 


2.69 


2.39 ± 0.16 


2.78 


DELPHI 


130.2 


9.7 ± 1.9 


8.1 


10.2 ± 3.1 


8.3 


136.2 


6.6 ± 1.6 


7.0 


8.8 ± 3.0 


7.2 


161.3 


3.6 ± 0.7 


4.5 


5.1 ± 1.2 


4.6 


172.1 


3.6 ± 0.7 


3.8 


4.5 ± 1.1 


3.9 


L3 


161.3 


4.59 ± 0.84 


4.4 


4.6 ± 1.1 


4.5 


172.1 


3.60 ± 0.75 


3.8 


4.3 ± 1.1 


3.9 


OPAL 


130.25 


9.0 ± 1.8 


8.0 


6.8 ± 2.0 


8.0 


136.22 


11.4 ± 2.1 


7.0 


7.2 ± 2.1 


6.9 


161.34 


4.5 ± 0.7 


4.4 


6.2 ± 1.0 


4.4 


172.12 


3.6 ± 0.6 


3.8 


3.9 ± 0.8 


3.8 


183 


3.46 ± 0.29 


3.45 


3.31 ± 0.32 


3.45 


189 


3.13 ± 0.16 


3.21 


3.53 ± 0.21 


3.21 
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Table 5: The best estimate of the r? = ^/M|, 95%CL upper Umit r]+, and the 95%CL lower limit on MsIT'^''^. 



Data Set 


r/ = ^/M|(GeV-4) r?+ (GeV^) 


95%CL MsjT^I^ 


DY-CDF 
DY-DO 
DYcombined 


(0.00 ± 0.15) X 10-11 0.22 X lO-n 
(0.00 ± 0.12) X 10-11 0.15 X 10-11 
(0.00 ± 0.10) X 10-11 0.14 X 10-11 


819 GeV 
906 GeV 
925 GeV 


LEPII hadronic 
LEPII leptonic 
LEPII b,c 
LEPII 


(0.45 ii:i8) X 10-11 0.82 X 10-11 
(-0.34 ^oies) X 10-11 1.01 x lO-n 

(1.63 -\fi) X 10-11 4.23 X 10-11 
(-0.44 ±0:29) X 10-11 0.74 X 10-11 


591 GeV 
561 GeV 
392 GeV 
607 GeV 


HERA 


(-0.050 ^ai?) X 10-11 0.47 x lO-n 


679 GeV 


ALL 
ALL w/o HERA 


(-0.026 -o°o76) x IQ-ii 0.13 x IQ-n 
(-0.0055 iaio) X IQ-" 0.14 x IQ-n 


939 GeV 
925 GeV 
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